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ABSTRACT 

We present high spatial resolution (0.4", ~ 3.5 kpc) PdBI interferometric data on 
three ultra- luminous infrared galaxies (ULIRGs) at z ~ 2: two submillimetre galaxies 
and one submillimetre faint star forming radio galaxy. The three galaxies have been ro- 
bustly detected in CO rotational transitions, either ^^C0(J=4— )-3) or ^^C0(J=3— )-2), 
allowing their sizes and gas masses to be accurately constrained. These are the highest 
spatial resolution observations observed to date (by a factor of ~ 2) for intermediate- 
excitation CO emission in z ~ 2 ULIRGs. The galaxies appear extended over several 
resolution elements, having a mean radius of 3.7 kpc. High-resolution (0.3") combined 
MERLIN- VLA observations of their radio continua allow an analysis of the star for- 
mation behaviour of these galaxies, on comparable spatial scales to that of the CO 
observations. This 'matched beam' approach sheds light on the spatial distribution of 
both molecular gas and star formation, and we can therefore calculate accurate star 
formation rates and gas surface densities: this allows us to place the three systems in 
the context of a Kennicutt-Schmidt (KS)-style star formation law. We find a differ- 
ence in size between the CO and radio emission regions, and as such we suggest that 
using the spatial extent of the CO emission region to estimate the surface density of 
star formation may lead to error. This size difference also causes the star formation 
efficiencies within systems to vary by up to a factor of 5. We also find, with our new 
accurate sizes, that SMGs lie significantly above the KS relation, indicating that stars 
are formed more efficiently in these extreme systems than in other high-z star forming 
galaxies. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: ISM 



1 INTRODUCTION 
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Recent evidence has suggested that the star formation activ- 
ity of the universe peaks at redshifts between 1 and 3, and 
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as such an understanding of the systems involved in forming 
stars at these cosmological redshifts is of great importance. 
Of particular interest arc the luminous submillimetre galax- 
ies, or SMGs (Small ct al . 2002 j , which peak strongly at 
z ~ 2 ( Chapman et al.||2005[ ). First detected with the Sub- 



millimeter Comm on-User Bolometer Array (SCUBA, Hol- 
land et ar]|l999 1, these sub-mm selected (SssoMm ~ 3 mjy) 
systems have prodigious star formation rates an order of 
magnitude greater than optically-selected galaxies at com- 
parable redshifts ( Si 1000 yr~^), lar ge reservoirs of 
molecular gas (~ 10^" Mq; 



Greve et al. 



20051, and are 



thought to contribute a substantial fraction of the total cos- 
mic star formation rate density at their peak epoch (~ 50% 
at z ~ 2; |Blain et aL]|1999[ |Blain et aL][2b02i ). Being sub- 



mm selected, the bolometric emission from these systems is 
dominated by cold dust (Td ~ 30 K), and they are relatively 
faint at 24/im and in the X-ray, indicating that the contri- 
bution from an AGN component is generally negligible. 

Observations of CO structures in SMGs have gener- 



ally suggested relatively small spatial extents ( Tacconi et al 



20081, typically having half light radii on the order 2 kpc. 
This, combined with the massive reservoir of gas in the in- 
terstellar medium (ISM), leads to a very high density (up to 
two orders of magnitude denser than their optically-selected 
counterparts - Nesvadba et al.||2007 l which triggers the in- 
tense star bursting behaviour that is characteristic of the 
SMG population. Many authors have concluded that the 
SMGs are therefore 'scaled up' versions of ultraluminous in- 
frared galaxies (ULIRGs, galaxies with Lfir > lO^'^ Lq), in 
the local Universe, being a population of compact star form- 



ing galaxies driven by major mergers (Sanders & Mirabel 



1996) with CO sizes comparable to the radii of the entire 
galaxy ( [Tacconi et aL|1999)|Iono et al. 12007] ). In this model, 
the star formation is concentrated entirely in the dense cir- 
cumnuclear region. Recent work has drawn upon evidence 
from features in the mid-IR (for example PAH emission in 
the study of 'Menendcz-Dclmestre et al."2009', and near-IR 
colours in the work of Hainline et al.^2009| ), using measures 
of the extinction to infer that the star forming region in 
SMGs is extended on scales of > 2 kpc; proportionally far 
larger than anything seen in local ULIRGs. It remains to 
be seen, therefore, whether conventional star formation pre- 
scriptions applicable to galaxies in the low redshift universe 
can account for the extreme star forming behaviour seen in 
SMGs. 

In this paper we report IRAM Plateau de Bure Interfer- 
ometer (PdBI) observations of three high-redshift star form- 
ing sources. Combined MERLIN -fVLA radio continuum ob- 
servations are used to analyse the spatially resolved star for- 



mation behaviour (e.g. Chapman et al.||2004 


Muxlow et al. 


|2005||Biggs & Ivison||2008| |Casey et al.||2009 


1, in combina- 



tion with the spatially resolved information on the molecu- 
lar gas content provided by the PdBI, to attempt to anal- 
yse the connection between the gas content and star forma- 
tion behaviour in these extreme high-redshift systems. The 
high-resolution of both the PdBI and MERLIN observations 
allow the sizes of the CO and star forming regions to be 
ascertained accurately and independently, allowing for ac- 
curate measurements of both the SFR surface density and 
the molecular gas surface density. As such, we can place 
our sources into the context of the Kennicutt-Schmidt law. 



which has been proposed as a possible 'universal star forma- 
tion law' ( [Bouche et al.|2007[ ). 

In §2, we outline the observations behind this work, 
and in §2.3 we outline the techniques used to analyse our 
data. In §3 we individually discuss our galaxies, in terms of 
their CO properties, star formation behaviour, and kinemat- 
ical parameters. §4 contains further discussion of merger vs. 
disc scenarios and star formation efficiencies for our galaxies, 
and attempts to place them in the context of a Kennicutt- 
Schmidt star formation law, and we conclude in §5. 

Throughout this work we assume a concordance ACDM 
cosmology with h = 0.71, Qa = 0.72, and rijvi = 0.28. 



2 OBSERVATIONS AND ANALYSIS 

We aim to compare matched beam (~ 0.4") observations of 
mid-excitation CO lines and radio continua in a selection of 
2 ~ 2 ULIRGs, two in the Hubble Deep Field North and one 
in the Lockman hole, whose optical counterparts have been 
identified as HDF132 {z = 1.999), HDF254 (z = 1.996), 
and Lockman 38 (z — 1.523) respectively. In this work they 
are referred to as such hereafter. This project was under- 
taken with a long-term goal of including similarly high res- 
olution of the resolved sub-mm continuum in order to fully 
characterise the nature of star formation, gas content, and 
dust in these extreme systems. We have chosen three sources 
with good CO detections which have been well detected with 
MERLIN radio observations. In particular, the radio prop- 
erties of the three sources were chosen to be representa- 



tive of the SMG radio population as a whole (see Chapman 
|et al.|,2004) , which has ~ 1/3 being extended in the radio 
(HDF254), and ~ 2/3 being radio compact (HDF132 and 
Lockman 38), so the large CO structures discussed below 
do not result from a selection bias. 

The three sources also span the complete range of z ~ 2 
ULIRG spectral energy distribution (SED) types - Lockman 
38 suggests a unusually cold SED, with a dust temperature 
Td ^ 20K. HDF254 su ggests a hot SED, with a Td of 65 ± 
14K ( Casey et al.|2009 ): it is strong in the radio, weak in the 
sub-mm, and has significant 70/im emission - it is, in fact, 
one of the few z > 1 ULIRGs (< 5%) detected by Spitzer at 
70/im. HDF132 is intermediate between these two sources, 
with an estimated Td of 40K ( [Chapman et al.|2005| ). 



2.1 PdBI CO observations 

We observed CO rotation lines in the sub-mm faint star 
forming radio galaxy (SFRG) RGJ123711. 344-621331.0 (= 
HDF254), the sub-mm galaxy SMMJ123618. 33-^621550. 5 
(= HDF132), and the SMG SMMJ105307-h572431.4 {= 
Lockman 38) during February and March 2009, using be- 
tween 4 and 6 of the antennas (of a total of 6). All sources 
were observed with the 'long baseline' A configuration. 

HDF132 was observed in ^2CO(J=4-)-3) on 2009 Febru- 
ary 15 and February 19, with on-source integration times 
of 8.9 and 2.9 hours respectively. Four of the six antennas 
were used. The observed frequency of the '^^C0(J=4— >-3) line 
was 154.004 GHz, putting the CO-derived redshift at z — 
1.999 ± 0.001. HDF254 was observed in i^CO(J=4^3) on 
2009 February 16 with a total integration time of 6.2 hours, 
using five antennas. ^^CO(J=4— >-3) was observed at 154.122 
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GHz, setting the source at z = 1.996 ± 0.001. Lockman 38 
was observed with all six antennas on 2009 February 20 and 
2009 March 12, with on-source integration times of 5.6 hours 
and 3.2 hours respectively. ^^CO(J=3— >2) was observed at 
154.122 GHz, setting the source at z = 1.523 ± 0.002. The 
synthesised beam size at 154 GHz is ~ 0.4" (see panel insets), 
~3.4 kpc at the redshifts of our sources. 

The SMG SMMJ123711.98+621325.7 {= HDF255) was 
also within the field, being just 7" south-east of HDF254. It 
was, however, undetected with our long baseline configura- 
tion observations, which can be attributed to a combination 
of primary beam attenuation effects, and the CO source be- 
ing sufficiently extended that all the flux is resolved out, con- 
sistent with a large radio size (~ 3.2 kpc) from MERLIN as 



reported by Chapman et al. (20041. We include HDF255 in 



our compilation of SMGs from the literature, however, and 
use the fact that it was undetected in our high resolution 
observations to place a lower limit on the source diameter 
of 1" (hence, all surface densities derived for HDF255 will 
strictly be upper limits). 

All sources were observed in good meteorological con- 
ditions, and typical system temperatures were 80-250K. 
Antenna efficiencies were typically 25-35 Jy K"\ Bright 
quasars (primarily 3C84 and 3C273) were used for fiux cali- 
bration. All data were reduced using the CLIC and mapping 
routines as part of the IRAM GILDAS package, and the re- 
sultant data cubes were rebinned to a frequency resolution 
of 20 MHz (a velocity spacing of 38 km s~^), before analysis 
with our own IDL routines. 



2.2 MERLIN 1.4 GHz continuum observations 

High resolution observations of the galaxies' radio continua 
were obtained using the Multi-Element Radio Linked Inter- 
ferometer Network (MERLIN; Thomasson|1986 1, combined 
with the VLA. Radio data for the HDF sources were ob- 
tained by 'Morrison et al. (prep), and data for Lockman 38 



were obtained by Ibar et al.| (2009) . 

The FWHM synthesised beam sizes for the maps are 
(semi-major axis vs. semi- minor axis in arcsec) 0.204 x 0.193, 
0.205 X 0.192, and 0.5 x 0.5 for HDF132, HDF254, and Lock- 
man 38 respectively. These beam sizes are all comparable in 
resolution to the long-baseline configuration beam size of 
the PdBI. The high spatial resolution of these observations 
allows for the star formation regions to be mapped in de- 
tail (rather than returning a single integrated value for the 
SFR), which is useful both for distinguishing point source 
AGN emission from diffuse star formation, and for compar- 
ing the extent and shape of the star forming region with the 
CO morphology (see Biggs fc Ivison||2008 1 . All three of our 
galaxies were well detected and fully resolved. 



2.3 Analysis 

Figs[l]|2] and|3]show the CO observations and spectra, com- 
bined with optical images and MERLIN -I- VLA radio contin- 
uum observations (as detailed in the captions). Tables[l]and 
[2]contain the observational results and other derived param- 
eters. All lines were successfully detected at a very high level 
of significance (> 6a). 

The two HDF sources (132 and 254) have optical images 



comprised of b,v,i &i z band images from the ACS, stacked 
with equal weighting. In the case of HDF132, NICMOS 1.6 
micron photometry (using the F160W filter) has been added 
for clarification. There is no ACS image for Lockman 38; in- 
stead, we use an I band image, taken with the CFH12K in- 
strument on the Canada- France-Hawaii Telescope (CFHT). 
Being ground-based, this is necessarily a lower-resolution 
image than the HST-ACS photometry, with a resolution of 
~ 0.7". 

The middle panels of Figs[l][2j and |3] also show the CO 
data colour coded by velocity. These were created as follows: 
each data cube has a channel range across which the signal 
to noise of the source is maximised. For each source, this 
channel range was then split into three equal bins, and the 
CO flux for each point on the map was integrated for each 
bin separately. This results in three separate maps of the 
source, corresponding to the three bins in channel space (or 
three bins in velocity space). These three maps were then 
colour-coded appropriately and stacked. 

The star formation rate given in Table [2] has been de- 
rived using the empirically ca librated FIR-radio correlation, 
( [2OO1I P" 



taken from 



Yun et al. 



Lfir = 47rD|, (8.4 x 10^*) S'1.4 (1 + ^ 



(1) 



where Dl is the luminosity distance in metres, 5*1.4 is the 
radio fiux density at 1.4 GHz (in W m^^ Hz~^), and a is 
the synchrotron slope used to K-correct the 1.4 GHz obser- 
vations to the appropriate source redshift (a is taken here to 
be 0.8). The Lfir-SFR conversion was taken from Kennicutt 
(1998a) : 

(2) 



SFR = 4.5 X 10"*"* Lfir (ergss"') 



It is important to note that several assumptions underlie the 
direct conversion between Lfir and a SFR; namely that all 
the UV radiation produced by young stars is absorbed by 
dust and re-radiated in the IR, and that the dust heating is 
brought about by a young stellar population (^30 Myr). If 
the stellar population is more evolved, the true SFR can be 
lower than that derived using the equation above by up to 
a factor of 2. 

CO fiux densities were estimated by fitting a Gaussian 
profile to the spectrum. From this, we calculate the CO lu- 



minosities using the relation given by Solomon et al. ( 19971 
L^o = 3.25 X 10' Sco (1 + z)-^ dI 



(3) 



(where 5*00 is the CO flux in Jy km s~^, i^obs is the fre- 
quency of the observed CO line in GHz, and the luminosity 
distance Dl is in Mpc), from which molecular gas masses 
were derived using a conversion factor appr opriate to local 
ULIR Gs of X = 0.8 Mq K km s'^pc^ (see 



2008 for discussion). 



Tacconi et al 



These molecular gas masses have been derived from CO 
luminosities resulting from moderate-excitation transitions 
(either 4— >-3 or 3— >-2); the extrapolation of an accurate gas 
mass, therefore, requires a derivation of a CO(l— >-0) flux 



^ The radio flux densities used here were obtained from the VLA- 
only data which recover the total flux density of each source, and 
are more sensitive to diffuse, extended regions of emission which 
may be resolved out by the longer baseline MERLIN data. 
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density. We use the brightness temperature conversions de- 



rived from SED modelling of SMGs by Weiss et al. (20071, 



Casey 



which leads to the following line ratios (as reported by 
et al.|2009| ): Sco3-2/Scoi-o = 7.0 ±0.5; Sco4-3/ScoTI^r^ 
10.0 ± 0.8. Most ULIRGs are thermalised up to at least 
the CO(3— >-2) transition, so using these mid-excitation lines 
(Juppcr = 3 or 4) should not significantly reduce the extent 
of the CO emission region when compared to observations 
of CO(l— >-0). Any size difference will be due to the fact that 
the higher order transitions trace the actively star forming 
gas, (which is of interest to us) whereas the lower transitions 
trace the cold gas either in the outer halo, or accreting onto 
the galaxy. 

Use of the long baseline 'A' configuration will inevitably 
resolve out flux when compared to observations made with 
more compact conflgurations. For HDF132 and HDF254 we 
compare our gas masses derived from our observations to 
those derived using the PdBI 'D' conflguration, to get an 
estimate of how much flux has been lost, which is typically a 
factor of ~ 2. For Lockman 38, which has not been observed 
in 'D' conflguration, we estimate a true gas mass based on 
the mean flux lost for the other two sources. We use the 
(more complete) 'D' configuration-derived gas masses when 
calculating gas surface densities for all systems. 



3 RESULTS 
3.1 HDF132 

HDF132 is both the brightest of our three sources at mil- 
limetre wavelengths, and the most highly obscured (being 
virtually undetected in the optical), implying an extremely 
high column density of dust. However, it is detected in the IR 
(as evidenced by the NICMOS F160W imaging, which was 



obtained by Conselice et al. 2008 as part of PID #11082), 



and is sub-mm bright, with a SCUBA 850/xm fiux of 7.3±1.1 
mjy. The strong 1.4 GHz continuum luminosity implies a 
high rate of extinction-corrected star formation. This is cal- 
culated (from the radio continuum, using the methods de- 
scribed above) to be 1500±100 Mq yr ~^ (here and through- 
out, quoted errors on the SFR represent fractional errors in 
the measured radio flux density, and do not take into ac- 
count other uncertainties, for example scatter in the LpiR- 
radio correlation). The velocity-coded image in the centre of 
Fig.[T]shows what appears to be either a disc-type structure 
or a late stage merger. 

Spectra taken at separate points along the source (see 
positional arrows in Fig. [T]) reveal the CO spectrum to be 
double peaked at all points, with the relative intensity of the 
peaks varying along the axis of the source. (The spectra have 
been more coarsely binned, to 40 MHz [~ 76 km s~^] for clar- 
ity). This could possibly be indicative of a 'beam smearing' 
effect, whereby the size of the beam (see inset ellipse) causes 
both velocity components to be detected. However, the posi- 
tion in velocity space of each peak is approximately constant 
across the galaxy, with just their respective strengths vary- 
ing. This, coupled with the relatively small beam size (0.4" 
- smaller than the source), suggests that the two peaks may 
be components of line of sight rotation. 

The position-velocity diagram (Fig. [l] bottom left) 
shows evidence for ordered gas motion, with two distinct 



clumps separated in velocity space; while they do overlap 
in physical space they are offset somewhat (resulting in the 
two peaks changing in intensity along the axis), this could 
be a signature of rotation (e.g. ,Funcs ct al. 2002) , but also 
may be indicative of a late stage merger. 

Fig. [5] shows the kinematics of HDF132 in more de- 
tail. Data were taken along the major axis of the source at 
~ 0.15" intervals (or ~ 1/3- 1/2 of the beam FWHM). The 
bottom panel shows the velocity at each point, obtained by 
fitting Gaussian curves to each peak of the source. To reach 
the final mean value for the velocity spread the positions of 
the two peaks were weighed by their respective strengths. 
It is likely that the value of the velocity shear obtained in 
this manner underestimates the true value, as beam smear- 
ing effects cause there to be significant emission from both 
components at most points. A better estimate for the true 
Vd would be Av/2, where Av is the separation of the two 
peaks in velocity space. For HDF132, Av/2 = 288 ± 19 km 
s -\ 

The middle panel of Fig. |5] shows the velocity disper- 
sion at each point, taken to be the la value of the Gaus- 
sian fit to the peak. The 'red' and 'blue' peaks have been 
treated separately, and are coloured red and blue respec- 
tively. The top panel compares the CO and the 1.4 GHz 
flux density at each point. The small size of the source (D 
< 1") only allows for four points to be taken with enough 
S/N to allow robust spectra to be taken. This precludes 
the full elucidation of the velocity behaviour that would be 
able to characterise the source as either a rotating disc or 
a merger. However, using the velocity dispersion (Fig. [5] 
middle panel) it is possible to estimate the lower limit on 
Vd/a (see Genzel et al. [200^ 1, a simple way of characterising 
the system as either 'rotation-dominated' (vd/cr > 2 — 3) 
or 'dispersion-dominated' {vd/a < l|^ The mean veloc- 
ity dispersion is 113 ± 19 km s^^ and 157 ± 32 km s~^ 
for the red and blue components respectively, leading to a 
Vd{smi)/a — 2.03 ± 0.26. This is, to first order, compara- 
ble to the value claimed for other high redshift star forming 
discs (e.g. For ster Schreiber et al.||2006 
[Law ct al. 200;^^ 

Given the fact that it is not apparently dispersion- 
dominated, and that it is homogeneous, compact, and non- 
clumpy on scales < 0.5", we identify HDF132 as a disc-type 
structure, as opposed to an early-stage major merger (al- 
though of course a late-stage major-merger remnant is very 
difficult to distinguish from a true secular disc given the 
limitations of our observations). 

The top panel of Fig. [5] shows a spatial breakdown of 
the CO and 1.4 GHz distributions. The two trace each other 
closely, both being centrally peaked and concentrated most 
strongly over ~ 0.5", and reaching a maximum value at 
the same point. The radio continuum emission peaks more 
strongly, however, falling off more sharply than the CO emis- 
sion either side of the central peak. Based on its CO mor- 
phology, we assign a source size (based on the CO emission 
haff light radiuiQ of 0.39" by 0.24", = 3.3 x 2.0 kpc at 
z = 1.996. 



Genzel et al. 112006 



^ For reference, z ~ discs typically have (vd/cr) ~ 10 — 20 | |Dib| 
|et al.|[2006[ l, whereas high-z star forming discs have (vd/cr) of 
order unity t o a f ew ( Genzel et al.|2008| Stark et al.|2008 i. 
Note that |KennicuttN199Sbl used the full measured e 



extent 
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Table 1. CO properties for HDF132, HDF254, and Lock38. Note that HDF132 and LockSS have two entries, for values derived from 
each of their double peaks. *H2 masses derived from compact-configuration data, to avoid resolving out fiux. 



Source 


RA 


DEC 




Sco 


L' CO 


L' CO(l-O) 


M(H2)* 




[J2000] 


[J2000] 




[Jy km s^-"-] 


[lO^f K km s-i pc2] 


[10"^ K km s-i pc2] 


[lO'J] 


HDF132 


12:36:18.337 


62:15:50.45 


1.996 


0.77 ±0.12 


0.94 ± 0.14 


0.59 ±0.11 


7.2 ± 1.5 








2.001 


0.50 ±0.08 


0.65 ±0.09 


0.40 ± 0.07 




HDF254 


12:37:11.354 


62:13:30.90 


1.996 


0.63 ±0.09 


0.81 ± 0.11 


0.51 ±0.09 


8.3 ± 1.7 


Lockman 38 


10:53:07.100 


57:24:31.87 


1.523 


0.51 ± 0.07 


0.68 ± 0.49 


0.53 ±0.36 


8.0 ± 3.9 








1.519 


0.12 ± 0.03 


0.16 ±0.04 


0.12 ± 0.03 





Table 2. Other properties for HDF132, HDF254, and Lock38. Note that FIR luminosities and SFRs are derived from the 1.4 GHz radio 
fiux densities. Note also that the SFR for HDF254 has been corrected for AGN contamination, as described in §3.2| The 1.4 GHz radio 
fiuxes were obtained by Morrison ct al. (prep). 



Source 


Sgso 


Si. 4 


L(FIR) 


SFR(FIR) 


log SsFR 


CO size 


(radii) 


1.4GHz size (radii) 


Mdyn 




[mjy] 




[10^2 Lq ] 


[M0 yr-i] 


[Mq yr-i kpc-2] 


[a"x 


b"] 


[a"x b"] Mq] 


[IOIOMq] 


HDF132 


7.3 ± 1.1 


172 ± 8.4 


8.6 ±0.6 


1500 ± 100 


2.19 


0.39" X 


0.24" 


0.24" X 0.18" 


21.8 ± 1.8 


HDF254 


2.5 ± 1.2 


126 ±8.7 


6.5 ±0.8 


1050 ± 90 


1.76 


0.52" X 


0.16" 


0.41" X 0.14" 


9.6 ±0.5 


Lock38 


6.5 ± 1.9 


56.2 ± 22.4 


1.5 ±0.6 


250 ± 80 


1.13 


0.56" X 


0.31" 


0.29" X 0.25" 


6.2 ± 1.1 



The long-baseline 'A' configuration observations often 
resolve out CO flux wlien compared to a lower-resolution 
observation (such as one made in compact 'D' configura- 
tion with the PdBI). To ensure that we deal with accurate 
gas surface densities, we take the source size from the high- 
resolution observations as described above, and the total gas 
mass from low-resolution (compact configuration) observa- 
tions. The total gas mass for HDF132 is 7.2 x 10^ Mq. (from 
as-yet unpublished PdBI observations |^ This leads to a gas 
surface density (taking the total gas mass from a mean of 
the two component peaks) of Eh2 ~ 540 Mq pc~^. To com- 
pare, the star-forming region (R1/2 derived from the radio 
source extent) is 0.24" by 0.18". The star formation surface 
density based on the FIR-inferred SFR is therefore Esfr = 
150 Mq yr kpc'^ 

Given that we have identified HDF132 as a disc-type 
structure, we can use a simple disc rotation model to cal- 
culate its dynamical mass: given its radius of 3.52 kpc and 
peak-to-peak rotation velocity of 288 ± 19 km s ~^ , we calcu- 
late a dynamical mass of Mdyn sin'^ i = 5.5 ± 0.5 x 10^'' Mq 
(it is important to note that this is assuming the gas is 
undergoing Keplerian motion). Adopting a fiducial incli- 
nation angle of 30°, we find the dynamical mass to be 
Mdyn = 2.2 ±0.2 X 10" Mf^ 



3.2 HDF254 



IQ. 



HDF254 is the second brightest (in the radio) of the three 
sources surveyed. It is also the most extended of the sources. 



of the emission to calibrate the KS law, whereas we use a 'half 
light' radius as standard for interferometric observations. Any 
bias introduced by this dilference should be negligible. 
■* This suggests that the high-resolution observations have re- 
solved out a factor of 1.8 in flux. 



having > 0.5 mJy CO emission extended over nearly 1.5" 
(= 12.6 kpc at 2: = 1.999). It broadly consists of two 'knots' 
of CO emission, an 'northern' knot to the north-east, and 
a 'southern' knot to the south west. In addition, there is a 
slightly weaker CO emission component ~ 0.5" to the west 
of the southern knot, with a CO flux of 0.67 mJy - this is a 
strong MERLIN point source (with a peak 1.4 GHz bright- 
ness of ~36 ^Jy/beam), which suggests the presence of an 



AGN. However, the Ha/NII observations of [Swinbank et aL 



( 2004[ ) , and PAH features in the NIR spectra taken by Pope 
et al.' (2008) both show little or no evidence for an AGN. 



The lower CO knot is detected in the ACS optical image 
(Fig. [2j left panel) , and also robustly detected in 1.4 GHz 
continuum (Fig. [2] middle panel), implying a high rate of 
concomitant star formation. By contrast, the northern knot 
is undetected optically and has a very weak detection in 
the 1.4 GHz continuum (although radio continuum emission 
is detected above the 3a sensitivity limit of ~ 10 fiJy at all 
points along the source). As the velocity-coded image shows, 
the two knots are distinct structures in velocity space. Spec- 
tra taken across the source show it to be single peaked, with 
a broad velocity dispersion in the northern knot, and dou- 
ble peaked with lower velocity dispersions in the southern 
knot. The middle spectrum of Fig. |2] shows a 1-dimensional 
CO spectrum taken at the position of the MERLIN point 
source. This shows either continuum emission or a broad 
peak between -1000 km s~^ and -200 km s~^ relative to the 
systemic velocity, at a level of ~ 0.75 Jy km s~^. Integrat- 
ing the spectrum between these two velocities (where the 
detection is most robust) gives a signal to noise ratio of ~3. 

The SFR, calculated using integrated low-resolution 
VLA measurements of the radio continuum, is 1220 ± 100 
Mq yr~^. This is possibly somewhat overestimated, how- 
ever, as the low resolution of the VLA will include the con- 
tribution of the possible AGN to the radio fiux, artificially 
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Figure 1. HDF132. Top Left: Stacked ACS image with overlaid CO emission contours (yellow) and NICMOS F160W contours (green). 
The yellow contours shown in the left hand panel delineate the 97, 98, 99, 99.5 and 99.9th percentiles of the velocity integrated CO data 
(the resultant levels in mjy/beam are 0.43, 0.50, 0.74, 1.18, and 1.90). Bottom Left: Position-velocity diagram taken along the major 
axis of the galaxy (PA = 228°; defined by the position angle of the best fitting ellipse to the CO map), with size and velocity scales 
shown inset. Centre: Velocity-coded CO emission, with overlaid MERLIN 1.4GHz radio contours (white; levels in /tjy/beam are 10, 12, 
34, 40, and 65). The inset velocity shift legend represents the centre of either the 'red' or the 'blue' bin, respective to the centre of the 
'green' central bin. The CO has been plotted starting at a significance level of 2(t (in the integrated map); all points with detections 
< 2tr have been left black. The radio beam has semi- major and semi-minor axes of 0.204" X 0.193" respectively. The outer contour of 
the inset ellipse (bottom left) corresponds to the FWHM of the CO beam, which is ^ 0.4". Velocities are coded relative to the systemic 
velocity, and are shown inset. Right: ID spectra, binned to 40 MHz, taken at the positions indicated. The spectra have been fitted with 
double-peaked Gaussian profiles. 



boosting the apparent SFR by up to ~ 1/3 - if the western 
unresolved MERLIN radio knot is indeed due to an AGN, as 
expect ed from the strong hard X-ray emission. [Casey et al.| 
1 2009| report a corrected SFR of 1050 Mq yr \ which we 
take to be the true SFR hereafter. 

Fig. |6] shows a kinematic breakdown for HDF254. All 
parameters were derived as for HDF132 above. The observed 
velocity shear of 80 ± 28 km s~^ (although again this is a 
lower limit due to uncertainty in the source inclination) was 
obtained by fitting a single peaked Gaussian profile to the 
CO emission at points along the major axis of the source. 
As with HDF132, this method can drastically underestimate 
the true velocity shear across the system, and in the case 
of a merging pair may have little physical interpretation 
(aside from shedding light on the line of sight component 
of their predominantly transverse velocities). The velocity 
profile (Fig. [6] bottom panel) resembles a rotation curve, 
but it may be that the low spatial resolution (~ 4 kpc beam 
size) of our observations is inadequate to resolve the char- 
acteristic signature of a merging pair of bodies. The clumpy 
(non-disc-like) structure apparent in Fig. [5] and the fact that 



the CO morphology changes from being single- to double- 
peaked as the two separate 'knots' are examined, both hint 
at the source being composed of two distinct bodies which 
are undergoing the early stages of a major merger. 

The position velocity diagram (Fig. [2) bottom left) 
shows the gas motions to be both disordered and complex, 
which is highly suggestive of a merging system. The two 
spatially separated bodies are clearly visible, being sepa- 
rated along the position axis and distinct in velocity space. 
In particular, the southern component clearly consists of two 
kinematically separate features (manifesting as the double 
peaked spectrum), whereas the northern component is made 
up of one single velocity structure. 

The MERLIN-I-VLA radio observations show the star 
formation (measured using the FWHM of the 1.4 GHz emis- 
sion map) to be spatially extended over a region ~ 0.4 " (= 
3.36 kpc) in length. There are five distinct peaks, with a 
maximum separation of ~1.5". The star formation forms an 
extended structure, and four of the peaks positionally trace 
the CO emission. It is, however, much more highly concen- 
trated in the southern, or double peaked knot, as Fig. [2] 
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Figure 2. HDF254. Top Left: Stacked ACS image with overlaid CO emission contours (yellow). The yellow contours shown in the left 
hand panel delineate the 97, 98, 99, 99.5 and 99.9th percentiles of the velocity integrated CO data (the resultant levels in mjy/beam are 
0.40, 0.46, 0.56, 0.67, and 0.85). Bottom Left: Position-velocity diagram taken along the major axis of the source (PA = 237°; defined 
by the position angle of the best fitting ellipse to the CO map), with size and velocity scales shown inset. Centre: Velocity-coded CO 
emission, with overlaid MERLIN 1.4Ghz radio contours (white; levels in /ijy/beam are 16, 18, 21, 24, and 33). The inset velocity shift 
legend represents the centre of either the 'red' or the 'blue' bin, respective to the centre of the 'green' central bin. The CO has been 
plotted starting at a significance level of 2a (in the integrated map); all points with detections < 2a have been left black. The radio beam 
has semi-major and semi-minor axes of 0.205" X 0.192" respectively. The outer contour of the inset ellipse (bottom left) corresponds 
to the FWHM of the CO beam, which is ~ 0.4". Velocities are coded relative to the systemic velocity, and are shown inset. Right: ID 
spectra, binned to 40 MHz, taken at the positions indicated. The spectra have been fitted with either single- or double-peaked Gaussian 
profiles (the fit with the lowest used in each case). 



shows. The top panel of Fig. [6] shows this in more detail: 
while the southern component (right hand peak) is traced 
very closely by a peak in the radio continuum (with the small 
0.1 " ofTset attributable to astrometrical errors), the north- 
ern component (left hand peak) has very little concomitant 
star formation. 

HDF254 has a high velocity dispersion, with the mean 
value of a being 249± 18 km s~^. This is far higher than the 
likely value of Viot , inasmuch as a value of Viot can be placed 
on such a chaotic system. Examining only the southern knot 
(which due to its double peaked CO spectrum appears to 
have a rotational component), the value of Av/2 is 149 ± 12 
km s"'^, leading to a (vd sini/a) of 0.67±0.11, implying that 
the southern knot is dispersion dominated, and it likely to 
be undergoing a major merger. While the northern knot is 
single peaked, and thus has an undefined rotational velocity, 
its velocity dispersion is high enough, at 285 ± 25 km s~^, 
to say with some confidence that the northern component 
too is dispersion-dominated. 

Treating the two knots separately it is possible to esti- 
mate the mass of the system. The northern knot is a single 



peaked source, and as such using the peak-to-peak velocity 
shift as a measure of the Keplerian rotation is not viable. 
Instead, we adopt a velocity-dispersion approach, as in |Erb| 
et al.|(|2006b|, whereby 



Mdyn = 



CcrV 

G 



(4) 



with the factor C being 3.4 for a rotating disc at an average 



inclination angle, and 5 for a sphere (see e.g. Erb et al. 
|2006b| ). The mean velocity dispersion of the northern knot 
region is 264 ± 24.8 km s~^. Assuming a C value of 4 (the 
body is most likely a turbulent ellipsoid, somewhere between 
the assumptions of a disc and a sphere), the dynamical mass 
of the northern knot is Mdyn = 6.9 ± 0.4 x 10^° Mq. 

As explained above, the southern knot has a signifi- 
cant line of sight rotation component, so it is possible to 
model the component assuming Keplerian motion. The knot 
has a radius of 0.2",which combined with the rotational ve- 
locity above gives a dynamical mass for the lower knot of 
Mdyn sin^ i = 5.6 ± 0.3 x lO'' Mq, or, using the fiducial in- 
clination of 30° as above, Mdyn = 2.2 ± 0.2 x 10^" Mq. 
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Figure 3. Lockman 38. Top Left: CO emission contours (yellow), over CFHT-CFH12K /-band (greyscale). The yellow contours shown 
in the left hand panel delineate the 97, 98, 99, 99.5 and 99.9th percentiles of the velocity integrated CO data (the resultant levels in 
mjy/beam are 0.23, 0.26, 0.32, 0.40, and 0.49). Bottom Left: Position-velocity diagram taken along the major axis of the source (PA 
= 193°; defined by the position angle of the best fitting ellipse to the CO map), with size and velocity scales shown inset. Centre: 
Velocity-coded CO emission, with overlaid MERLIN 1.4Ghz radio contours (white; levels in /ijy/beam are 7, 8, 9, 10, and 13). The inset 
velocity shift legend represents the centre of either the 'red' or the 'blue' bin, respective to the centre of the 'green' central bin. The CO 
has been plotted starting at a significance level of 2a (in the integrated map); all points with detections < 2a have been left black. The 
radio beam has a radius of 0.5". The outer contour of the inset ellipse (bottom left) corresponds to the FWHM of the CO beam, which 
is ~ 0.5". Velocities are coded relative to the systemic velocity, and are shown inset. Right: ID spectra, binned to 40 MHz, taken at the 
positions indicated. The spectra have been fitted with double-peaked Gaussian profiles. 




Figure 4. HDF255. Left: Stacked ACS image with overlaid CO emission contours, taken with the compact 'D' configuration with PdBI 
(yellow). For reference, the source in the top right is HDF254. Right: CO spectra, also taken in compact 'D' configuration, at the centre 
of the source. For reference, the strongly detected source to the NW of HDF255 (which is the source at the centre of the image) is 
HDF254. 
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Figure 5. Kinematic analysis for HDF132 
(SMMJ123618.33+621550.5). Spectra were taken at points 
along the major axis of the galaxy with ~0.15" spacing (~ 1/3 
- 1/2 of the beam FWHM), and resultant kinematic parameters 
were derived. Bottom panel: circular velocity; as the CO spec- 
trum for HDF254 is double peaked, the velocity was obtained by 
taking the mean of the two peaks, with the weighting of each 
component taken as the height of the respective peak . Middle 
panel: velocity dispersion, obtained via the width of the same 
Gaussian peak. The solid blue line shows the blue shifted peak, 
whereas the dashed red line shows the red shifted peak. Top 
panel: CO surface brightness (blue solid line). For reference, the 
MERLIN-obtained 1.4GHz continuum flux density measured at 
the same points has been overlaid (black dashed line and right 
hand axis). The dashed horizontal line is drawn at lO^Jy/beam 
(~ 3 times the image rms).. 



For comparison, it is possible to fit a single Gaussian pro- 
file to the spectra and derive the dynamical mass using 
the velocity dispersion as above. Doing so (the mean ve- 
locity dispersion in the lower knot is 201 ± 15 km s~^), and 
adopting a 'disc' C value of 3.4, gives a dynamical mass of 
Mdyn = 3.1 ± 0.2 X 10^° Mq. This is comparable (within a 
factor of two) to the value derived using the peak-to-peak 
shift: the difference may be due to the poor fitting and ar- 
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Figure 6. Kinematic analysis for HDF254. Spectra were taken 
at points along the major axis of the galaxy with ^O.IS" spacing 
(~ 1/3 - 1/2 of the beam FWHM), and resultant kinematic pa- 
rameters were derived. Bottom panel: circular velocity, obtained 
using the centroid of the Gaussian profile fit to the CO emission 
line. Middle panel: velocity dispersion, obtained via the width of 
the same Gaussian peak. Top panel: CO surface brightness (blue 
solid line). For reference, the MERLIN-obtained 1.4GHz contin- 
uum flux density measured at the same points has been overlaid 
(black dashed line and right hand axis). The dashed horizontal 
line is drawn at lO^Jy/beam (~ 3 times the image rms).. 



tificial broadening of the single Gaussian profile fit to the 
double peaked spectrum. The total dynamical mass of the 
entire system, therefore, is Mdyn = 9.6 ±0.5 x 10^° Mq. The 
dynamical masses of the separate components clearly mark 
the system as being a major merger, with a mass ratio of 
(2.2 - 3.1):1. 

To calculate surface densities, we integrate values across 
the surface of the entire elongated body of the source. As for 
HDF132, we take the total gas mass from a low-resolution 
observation of the same source, to avoid resolving out flux. 



Casey et al. ( 2009 \ report the total gas mass of this source 
to be 8.3 X 10^ Mq, suggesting that our high-resolution ob- 
servations resolve out a factor of 2 in CO flux. Based on 
the CO morphology, we assign source semi-major and semi- 
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Figure 7. Kinematic analysis for Lockman 38. Spectra were 
taken at points along the major axis of the galaxy with ~0.15" 
spacing (~ 1/3 - 1/2 of the beam FWHM), and resultant kine- 
matic parameters were derived. Bottom panel: circular velocity, 
obtained using the weighted mean of the two Gaussian peaks fit 
to the CO emission line. Middle panel: velocity dispersion, ob- 
tained via the width of the same Gaussian peaks. The solid blue 
line shows the blue shifted peak, whereas the dashed red line 
shows the red shifted peak. Top panel: CO surface brightness 
(blue solid line). For reference, the MERLIN-obtained 1.4GHz 
continuum flux density measured at the same points has been 
overlaid (black dashed line and right hand axis). The dashed hor- 
izontal line is drawn at 11/jJy/beam (~ 3 times the image rms). 



minor axes (measured using the half light radius as explained 
above) of 0.52" x 0.16", which leads to a gas surface density 
value of Eh2 = 513 Mq pc~^. Adopting an identical method 
for the SF surface area, we fit the radio emission area with 
semi-major and -minor axes of 0.41" x 0.14". This gives a 
SFR surface density of Esfr — 58 Mq yr ~^ kpc~^. 

3.3 Lockman 38 

Lockman 38 is the weakest radio source in our sample, sug- 
gesting that it represents an example of the extreme radio- 



faint end of the SMG population, having a SFR typical of 
the faintest 20% of SMGs as presented by |Ivison et al.|2007| 
This is similar to the brightest BX/BzK galaxies at 2 ~ 1.5. 

The CO emission is moderately extended in the north- 
south direction, with the strongest emission coming from 
an elliptical region 1.1" in diameter, 9.3 kpc at z — 1.523. 
It is comprised of two main velocity components, a small 
blue-shifted region to the north and a larger (and brighter) 
red-shifted region to the south. As with HDF132, there is a 
fainter (by a factor of ~2) 'tail' emission component, with 
a velocity intermediate between the two 'main' regions. ID 
spectra taken at positions along the source show it to be 
double peaked at all points. As with HDF132, the position of 
each respective peak is close to constant, with the position of 
the red-shifted peak varying by ±33 km s~^ along the source, 
and the position of the 'blue' peak varying by ±38 km s~^. 
Unlike HDF132, however, the peaks are highly asymmetric, 
with the red-shifted peak near ?; ~ km s~^ being around a 
factor of 5 brighter than the blue-shifted peak at —400 km 
s-i. 

The position-velocity diagram shows the gas motion 
to be somewhat disturbed. While there is no evidence for 
spatially separated bodies, there are two distinct kinematic 
components, which are separated by u ~ —400 km s~^ 
in velocity space. These peaks are asymmetric (as dis- 
cussed above), and do not display the signature of a disc- 
component. The position-velocity diagram suggests a late- 
stage merger interpretation of Lockman 38 might be most 
appropriate, where the components have coalesced to form 
a single body, but are still kinematically distinct. 

The MERLIN±VLA radio continuum map shows the 
1.4 GHz emission to be strongly peaked within the bright- 
est region of the galaxy, with little or no detected emission 
in either the blue-shifted region to the north or the fainter 
tail. A spatial breakdown of this (Fig. [7j top panel) shows 
the radio emission falling below the sensitivity threshold of 
10/iJy for all but the inner ~ 0.3" of the source. 

A kinematic investigation of the source (Fig. [7| shows 
it to have the largest velocity gradient of the three sources 
(bottom panel), having a maximum Vrot, obtained as above 
using a weighted mean of both peaks, of 164 ± 24 km s~^. 
This velocity varies smoothly across the source. Using the 
peak-to-peak shift {Av/2), a better measure of the true max- 
imal rotation velocity, v^ot = 229 ± 28 km s~^. The velocity 
dispersion for the two separate peaks is shown in Fig. [7] The 
value of a for both peaks varies along the source, though in 
opposite directions. The mean value for the dispersion is 
79 ± 22 km s"^ and 123 ± 24 km s"^ for the blue-shifted 
and red-shifted peaks respectively. The value of Vd{sini)/a 
for this source (taking a mean of the two peak dispersions) 
is 2.3 ± 0.14, placing the source very much in the rotation- 
dominated regime. However, due to the disordered and non- 
rotational gas motions revealed in the position-velocity di- 
agram, we class this source as a late-stage merger, as it is 
not dispersion-dominated, and apparently non-clumpy on 
the scales revealed by our beam, yet does not seem to be- 
have kinematically as a true disc. 

As described in §2. 3| above, Lockman 38 has never been 
observed in the compact configuration required to measure 
a truly accurate gas mass - we can, however, estimate the 
true value by assuming that with our high-resolution obser- 
vations we resolve out a similar fraction to the other sources 
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in our programme. We therefore assume that we lose a fac- 
tor of 1.9 in flux, and adjust our derived gas surface density 
accordingly. 



Ibar et al. ( 2009 1 VLA imaging of Lockman 38 reveals 



a total flux density of 32.3 ± 7.2 /ijy, if the source size is 
assumed to be identical to that of the beam. However, the 
source is faintly resolved (vaguely consistent with a bipolar 
morphology) even by the 3.9" beam. If the source size is 
deconstrained, the total flux density is 56. 2± 22.4 /ijy, which 
is taken to be the true 1.4 GHz flux hereafter. This large 
radio extent (resolved by the large beam) suggests that in 
addition to the compact core there is a low-level extended 
starburst being fuelled by the extended CO reservoir. The 



610 MHz flux density, obtained by Ibar et al. ( 2009 1 with 



the GMRT is 59.8 ± 17.9 ^l^y, which supports the model of 
this system as a starburst (as opposed to a steep-spectrum 
AGN). 

Lockman 38 has an unusually low radio-derived star 
formation rate for an SMG, at 250 ± 80 Mq yr"^ This is 
despite having a sub-mm flux of 6.5 ± 1.9 mjy, implying a 
SFR of ~600 Mq yr-^ This low radio/sub -mm ratio is un- 
usual given its redshift, and hints at an unusually cold dust 
temperature for this source. Assuming the SFR(radio) to 
be correct, we again calculate gas and SF surface densities 
for the source. Due to the highly asymmetric nature of the 
peaks, we do not take a 'mean' value of L'(CO), but instead 
use the more significantly detected stronger peak. We fit the 
significant (R1/2) CO emission with an ellipse, with semi- 
major and semi-minor axes of 0.31" and 0.56"respectively. 
This gives a gas surface density (including the estimated fac- 
tor of 1.9 to account for flux resolved out) of Eh2 = 300 Mq 
pc~^. The radio continuum emission is more highly concen- 
trated than the CO emission, and is well fit by a 0.29" x 
0.25" ellipse; using the SFR given in Table [21 we arrive at a 
SFR surface density of Esfr = 14 Mq yr kpc~^. 

To estimate the dynamical mass of Lockman 38, we 
assume rotational motion, and use the peak-to-peak velocity 
shift as a measure of its maximal rotation velocity; we take 
the radius to be 0.56" (= 4.75 kpc) and the rotation velocity 
as 229±28 km s~^. Assuming Keplerian motion, we calculate 
the dynamical mass to be Mdyn sin i = 1.6±0.9x 10^" Mq. 
Taking the fiducial inclination of 30° , the dynamical mass is 



calculated to be M, 



dyn 



1.2 ± 1.1 X 10'" Mf7 



0- 



3.4 HDF255 

HDF255 falls within the field of view for the HDF254 ob- 
servation, being just ~7" to the south east. It is, however, 
undetected by our long-baseline 'A' configuration observa- 
tions, suffering from both resolution-based flux loss in the 
extended source, and primary beam attenuation effects (be- 
ing off phase centre). FigureWlshows ^■^C0(J=4— >-3) observa- 
tions taken with the low-resolution, compact 'D' conflgura- 
tion of the PdBI (obtained for our own pilot observations for 
this programme), along with stacked b,v,i &i z band optical 
imaging from the ACS. 

HDF255 has been observed in the radio continuum 
with both MERLIN and the VLA, so calculating the star 
formation rate surface density is possible. The FWHM of 
the MERLIN-I-VLA image has an elliptical size of 0.56" 
X 0.32", and the 1.4 GHz flux density observed from the 
VLA is 53.9 ± 8.1iiJy. K-correcting and applying the radio- 



far IR correlation (as per Eq. 1), this gives a SFR of 
621 ± 93 Mq yr~^, with a resultant SFR surface density of 
EsPR = 16 MQyr-^kpc-^ 

The galaxy is well detected in CO (~ 4.5(t) by the com- 
pact conflguration, but is poorly detected by the long base- 
line observations, suggesting that it is a low surface density 
source with a large extent which has been driven below the 
sensitivity threshold by the small beam (0.4") of the 'A' 
conflguration. The compact-conflguration flux for HDF255 
is Sco4-3 ~ 0.53 ±0.11 Jy km/s, which leads to a gas mass 
of 3.5±0.7x 10^ Mq (after correction to an L'co(i-o) equiv- 
alent as above). This relatively low gas mass, coupled with 
a source size larger than the 'A'-configuration beam at 154 
GHz, leads to the source being resolved out (we would only 
detect this were all the fiux to be concentrated within a beam 
size). We assume an approximate (lower limit) CO source ra- 
dius of 0.6", based on the high-resolution radio source size 
and the extent of the beam. This is a rough approximation 
at best, however, and should be treated with caution. This 
leads to a value of Egas ^ 125 Mq pc~^. 



4 DISCUSSION 

4.1 Mergers vs. discs 

It is challenging to distinguish between disc-type structures 



and late stage merger remnants at high redshift (e.g. Dasyra 
et al.|2008| [Melbourne et al.|2008| ). Indeed, SMGs are gener- 
ally classed as either undergoing a major merger, or consist- 



ing of a compact rotating merger remnant ( Swinbank et al 



2006[ [Tacconi et al.[[2008| [Shapiro et al][2068[ ). In our sam- 
ple, certainly HDF254 is a system undergoing an early stage 
major merger: as explained above in §3.2| the galaxy appears 
clumpy on sub-arcsecond scales, and has an extremely high 
velocity dispersion (see Fig. [6|, indicative of the chaotic and 
disordered gas motions created by such a merger. 

HDF132, however, does not appear to be undergoing an 
early stage merging event; its symmetrical velocity peaks, 
low velocity dispersion (relative to the rotational velocity), 
and ordered gas motions lend some credence to the non- 
merging picture. The NICMOS IR morphology of HDF132 
reveals an elliptical structure, with the major axis coincident 
(within the astrometrical errors) with the major axis of the 
CO emission. If we are to class HDF132 (and potentially 
Lockman 38) as a system not undergoing a major merger, 
or consisting of a recent merger remnant, we must be able 
to place these systems into a valid, cosmologically-motivated 
paradigm of galaxy formation which allows for extreme star 
formation rates and extended disc-like structures to emerge 
from hot /cold flow accretion type models. 

Recent simulations of giant galaxies at high redshifts of 
2 ( [Dave et al.][2009[ [Dekel et al^[2009[ ) have suggested 



that extended CO discs can form entirely secularly, without 
a major-merging event. In this 'smooth accretion' model, the 
systems formed are 'harassed' massive galaxies, which are 
being rapidly formed by filamentary accretion of material 
from the IGM. They can be pictured as undergoing a series 
of minor mergers (with dark matter halo ratios of >1:10), 
which are observationally indistinguishable from a smooth- 
accretion driven secular evolution. 



Narayanan et al. ( 2009 1 carried out simulations of the 
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Table 3. Properties of all SMGs, both observed in this programme and archival. All parameters for archival sources have been re- 
calibrated with the same conversion factors as given in this work. Half light radii for our sources are approximations to our (non-circular) 
sources for comparison only - see Table 2 for accurate source dimensions, which were used to derive all parameter surface densities. 



Source 


z 


Ri/2 (CO) 


Ri/2 (1.4 GHz) 


log Sgas 


log SspR MERLIN data reference 






[kpc] 


[kpc] 


[M0 pc-2] 


[Mq yr-i kpc-2] 








Our sources: 


















HDF132 


1.99 


2.7 


1.8 


2.73 


2.16 








HDF254 


1.99 


2.9 


2.3 


2.71 


1.76 








Lockman 38 


1.52 


3.7 


2.3 


2.48 


0.76 








Archival sources: 


















SMMJ105141-I-571952 


1.21 


3.0 


2.7 


2.59 


1.45 


Unpublished 


SMMJ123634-I-621241 
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star forming molecular gas within SMGs, using the a pnori 
assumption that SMGs are merger-driven systems. The ve- 
locity field (Fig. 6 of their paper) shows, in general, a chaotic 
and disordered structure at all times apart from a short pe- 
riod (~ 0.03 Gyr) in the late inspiral phase of the merger; a 
system at this stage is yet to undergo the 'final coalescence' 
into a compact merger remnant. It is this final coalescence 
that boosts the star formation rate and 850/im flux, which 
causes the system to be classed as an SMG, or more gen- 
erally a ULIRG. The reasonably bright 850/xm flux (> 7 
mjy) and high star formation rate (> 10^ Mq yr~^) argue 
against HDF132 being an inspiralling merger, suggesting in- 
stead that the flnal coalescence has already occurred. This 
conflict lends weight - at least from a simulation-led perspec- 
tive - to the argument that the behaviour of HDF132 cannot 
be easily explained by a merger. It must be noted, however, 
that while these simulations do manage to reproduce the 
high star formation rates seen in SMGs (> 10^ Mq yr~^) 
from a hydrodynamical perspective, the explicit model em- 
ployed is unable to account for the features of the entire 
SMG population, and as such must be treated with care 
when dealing with individual cases. 

4.2 Size of CO emission regions 

Recent work investigating the sizes of massive high redshift 
galaxies has shown that some high-z galaxies do, in fact, 
possess extended morphologies. For example, |Mancini et al.| 
( 2009[ ) investigated a sample of massive early type galaxies, 
flnding that their optical radii were comparable to the radii 
of similar mass galaxies in the z = Universe, suggesting 
that the low S/N of massive high-z galaxies causes their 
half-light radii to be systematically underestimated. 

HDF254, the only unequivocally merging system in our 
study, is comprised of two or three smaller bodies, each ~ 2 
kpc in diameter, comparable in extent to the mean size of 



SMGs identified by (for example) Tacconi et al. (20081. This 



is the first size measurement for a sub-mm faint ULIRG. 
However, the total system size is > 12kpc, and unlike other 



comparable mergers from Tacconi et al. ( 2006 1 and 



Tac- 



coni et al. (20081 such as HDF242, and 0931, there is a 



extents larger than typically seen in high resolution studies 
of SMGs. HDF132 is the most compact of our sources, but 
with a diameter of ~ 0.7" (= 5.9 kpc) is more extended 
than all but one of the 11 sources presented by |Tacconi| 



et al. 



( 2006 1. Lockman 38 has a diameter of 1.1" (9.5 kpc at 



= L523). The prevailing paradigm of SMGs being major 
mergers leads naturally to the conclusion that their molec- 
ular gas will be concentrated in a compact nuclear region: 
tidal torquing induced by the merger drives cold molecular 
gas into the centre of the galaxy, where it is heated, becomes 
gravitationally unstable, and triggers a starbursting event. 
This suggests that while systems undergoing major mergers 
may comprise a part of the SMG population, it is perhaps 
not exclusively so. 

4.3 The L'co-LpiR Relation and Star Formation 
Efficiency 

A common way to conceptualise the relationship between 
the gas content and star formation rate of a galaxy is in 
terms of a star formation efficiency, defined as the SFR per 
unit of star forming gas - i.e. SFR/M(H2). Measuring this for 
high- 2: galaxies requires an assumption of a CO-H2 conver- 
sion factor; while this is a necessary evil for calculating gas 
masses and gas surface densities, it is possible to frame the 
star formation efficiency in terms of observables (assuming 
the FIR-radio correlation), and use the continuum-to-line 
luminosity ratio Lfir/L'co- 

We have examined this continuum-to-line luminosity ra- 
tio for three 'sub-regions' of each of our sources, to better 
show the star forming behaviour within each. Fig. |8] shows 
the result for each of the galaxies. For comparison, other 
global values of Lfir/L'co for SMGs in the literature have 
been plotted - see the caption for references. All the sub- 
clumps scatter (by a factor of ~ 5) around the canonical 
Lfir/L'co relation as defined by Greve et al. ( 2005 1. This re- 



lation has a slope less than unity, indicating that as Lfir in- 
creases, as does the Lfir/L'co ratio: star formation is more 
efficient in ultraluminous systems like SMGs than in LIRGS 
and normal star forming galaxies. However, |Krumholz fc| 
McKee ( 2005 1 have suggested that it is not a star formation 



'bridge' of CO emission connecting the interacting clumps. 
Both HDF132 and Lockman 38, however, have CO spatial 



efficiency issue, but rather results from a disconnect between 
the gas being traced and the gas actively involved in star for- 
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Figure 8. CO luminosity plotted against FIR luminosity for 'sub- 
clumps' of our three sources. For reference, results for other SMGs 



in the literature have been plotted ||Genzel et al. 2003 Neri et al. 


2003 


Downes & Solomon"2003| |Sheth et al.| 


2004, iGreve et al. 


2005 


Daddi et al.||2009, Knudsen et al.||2009 


1. The dashed line 



represents the fit to SMGs and ULIRGs from Greve et al. (2005) 
of (log L'co = 0.62(log Lpipi) + 2.33). All our sources scatter 
around the canonical relation by up to a factor of 5. 



mation (see |Gao et al.|[2007[ for example). In this picture, 
CO emission with Jupper 5; 4 traces the ambient gas den- 
sity (which is non-hnear with SFR), whereas emission from 
molecules that trace the densest star forming cores (such as 
HCN) is better suited for measuring the density of the gas 
actually involved in star formation. Indeed, in the local Uni- 
verse L(HCN) is seen to increase linearly with L(FIR). The 
sub-clumps of our three SMGs follow this (non-linear) trend, 
and have Lfir/L'co ratios consistent with other SMGs. As 
can be seen from Fig. [8] Lockman 38 is the least-efficient 
SMG studied to date, further suggesting that it is a highly 
unusual system. 



4.4 Gas fractions and the uncertain CO 
conversion factor 

The question of the appropriate CO/H2 conversion factor 
to use in extreme high-z systems is an ongoing one. The 
generally accepted values for local ULIRGs are within the 



range 0.8 — 1.6 Mq K km s~ pc , substantially lower than 
the calibration used for 'normal' systems - this value (nor- 
mally referred to as the Milky Way, or galactic, conversion 
factor) is generally taken to be 4 — 5 Mq K km s~^pc^, and 
is derived both theoretically and empirically (see Appendix 
A of lTacconi et al. 2008 and references therein for discus- 
sion). Obtaining the correct value for any given system is 
difficult, however, as the true value is a complex function 
of many physical parameters, including the gas metallically 
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Figure 9. Dynamical mass plotted against molecular gas mass 
for all sources presented in this work (filled red circles). For com- 
parison, high-2 star forming galaxies and SMGs from Tacconi et 
al. (2008) and Greve et al. (2005) are also shown. Open red circles 
represent the position of the three sources presented in this work, 
if a Milky Way CO conversion factor is used to calculate the gas 
mass. Dashed diagonal lines, from bottom to top, represent con- 
stant gas fractions (i.e. Mgas/Mdynamical) of 0.1, 0.3, 0.5 and 0.7 
respectively. 



(Israel 20051, and the column density and temperature of 



the in tersteller gas ( |Sakamoto||1996| |Weifi et al.]|2001[ [Wall 
2007| ). 

The theoretical justification for using a significantly 
lower conversion factor for ULIRGs than normal star- 
forming galaxies hinges on the nature of merger-induced 
star formation. Whereas 'normal' star formation occurs in 
gravitationally bound molecular clouds (which exist as well 
defined separate structures in the ISM), the torque in- 
duced by a merger event creates a more centralised, homo- 
geneous molecular ISM, in which pressure plays a signifi- 
cant role in triggering a star forming event. This difference 
(gravitationally-confined vs. pressure confined star forming 
gas) theoretically justifies the use of a different conversion 
factor. Empirically, it is found that using a standard galac- 
tic conversion factor for ULIRGs often results in derived gas 
masses greater than their dynamical masses, necessitating 
the adoption of a lower value. 

The paradigm of SMGs being 'scaled up' versions of 
local ULIRGs has led to these lower values of Xco being 
adopted by most authors (i.e. X = 0.8 by Greve et al.|2005 
X = 1 by [Tacconi et al.|[2008| , which has led to a pre- 
dominance of low assumed gas fractions for SMGs, which 
typically have values of fgas (= Mgas/Mdynamicai) between 
0.1 and 0.3. Fig. [9] shows the relationship between dynam- 
ical mass and gas mass for the sources in this work, along 
with a selection of SMGs from the literature - diagonal lines 
of constant gas fraction have been drawn for clarity. A ma- 
jority of the galaxies shown (9/17) have low gas fractions 
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of < 20%, a value which would be modest for local star 
forming discs (e.g. Karachentsev et al.|2004 Bothwell et al. 
[2009) and seems utterly incongruous with the exceptional 
star formation rates achieved by SMGs. 

For the 3 ULIRGs considered in this paper, we have 
measured relatively extended CO structures of low gas sur- 
face density, comparable to local starburst galaxies, and 
distinctly different from the dense, nucleated gas of local 
ULIRGs. It is therefore reasonable to consider that the use 
of the CO conversion factor appropriate for local ULIRGs 
{X ~ 1) might not be justified in the types of systems un- 
covered in this work. We also show in Fig. [o] (as open red 
circles) the effect of using a Milky Way conversion factor 
{X = 4.3 Mq K km s~^pc^) for our three galaxies. 

A full theoretical and empirical treatment of the CO/H2 
conversion factor lies beyond the scope of this work. We 
do, however, include in the following discussion of the star 
formation law the possibility that the assumed value of X = 
0.8 Mq K km s~^pc^ may be significantly underestimating 
the true gas mass. 



-3 
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4.5 The Star Formation Law- 
There have been many recent studies on star formation at 
high redshifts, putting the extreme systems seen at z ~ 2 
(where the cosmic star formation rate density peaks) into 
the context of a 'universal star formation law'. The most 
common prescription, at both low and high redshifts, is the 
power-law scaling between the surface density of gas (Egas) 
and star formation (Espr) known as the Kennicutt-Schmidt 
(KS) law ( [Schmidt 1959; K ennicutt|1998b[ ). 

A full and correct calibration of the high-z KS law is 
crucial in order to facilitate studies of star forming galaxies 
at cosmological redshifts. With current instrumentation it 
is only possible to measure CO masses for the most lumi- 
nous systems: this has the result that all galaxies currently 
observed in CO aX z > 2 are ULIRGs, which are clearly 
unrepresentative of the high- 2: population as a whole. In ad- 
dition, obtaining accurate gas masses via CO measurements 
is highly time consuming, requiring on-source integration 
times of several hours; as such, obtaining well constrained 
gas masses for a statistically significant sample of high-z star 
forming galaxies would be prohibitively expensive in terms 
of instrument time. One way around this is to infer a gas 
mass based on the (more easily measurable) Ha emission 
size and star formation rate (e.g. [Erb et aL][2006a| . This 
requires the assumption of a universal star formation law, 
which has only been truly calibrated at low redshifts, and 
may be inaccurate for both the 'normal' (L* galaxies) and 
extreme systems (SMGs, SFRGs etc.) lying at cosmologi- 
cal redshifts. Clearly, for this kind of large scale work to be 
undertaken there is a need for high redshift galaxies to be 
placed into the context of a KS-type star formation law. 



Erb et al. (2006a I established a plausible KS-law rela- 



tion for BX/BM galaxies by inverting the KS law for gas 



surface densities. Tacconi et al. (20061 explicitly measured 



the gas surface densities in SMGs, and inferred the basic KS- 
law properties of SMGs, while [Bouche et ar] ( [2007[ ) reviewed 
the literature on this topic, and placed all these galaxies 
into a self-consistent framework, re-iterating the findings 
of both these authors. [Bouche et al.[ ( [2007[ ) report that, in 
general, extreme star forming galaxies at high redshifts lie 



log i:„j5 [M_sun pc ] 

Figure 10. Log gas surface density plotted against log star forma- 
tion surface density (the 'Kennicutt-Schmidt relation'; see |Ken^ 
|nicutt|[l998b}. Black crosses and green triangles represent the 
normal spirals and circumnuclear starbursts (respectively) pub- 
lished by Kennicutt (1998b). Blue squares denote previously de- 
rived values of Sgas and Ssfr for SMGs in the literature, and 
red circles denote the same parameters for the three new galax- 
ies presented in this work. Previously published SMG results 
have been adjusted to a common CO/H2 conversion factor of 
0.8 Mq K km s^'^pc^, and all SMG molecular gas masses have 
had a helium correction factor of 1.36 applied. Horizontal 'error 
bars' on the SMGs/SFRGs extending to the right show the range 
of gas surface densities achievable by varying the CO/H2 con- 
version factor, up to a Milky Way value of 4.3. The dashed line 
is the original {N = 1.4) slope given by Kennicutt (1998b), and 
the solid fit is the improved fit (A'^ = 1.71) applicable to high-z 
galaxies as presented by Bouche et al. (2007). 



above the 2: ~ KS law originally proposed by |Kenniciitt] 
(1998b), a fact which they attribute to an environmental 
enhancement of the star formation rate most likely due to 
collisions or pressure. When including the high-z galaxies, 
they modify the canonical KS slope of A'^ ~ 1.4, finding 
that A'' = 1.71 ± 0.05 provides a better fit. Even taking 
this steeper relation into account, the SMGs in their study 
still lie above the relation, indicating the very high efficiency 
of their star formation. [Bouche et al. ( 2007 1 were, however, 
forced to make assumptions, including the fact that the sizes 
of the CO region and the star forming region were identi- 
cal, using the half-light radius (R1/2) of the CO emission for 
both. 

Our high-resolution, matched beam study of SMGs al- 
lows us to improve over previous studies by taking indepen- 
dent observations of both the sizes of the star forming and 
the CO emission regions, and by utilising accurate and in- 
dependent measures of star formation rates and gas masses. 
As discussed in §3 above (and see Table [2|, we find that 
in all instances the star forming region is centrally peaked, 
and is extended over a smaller region than the CO morphol- 
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ogy. Taking the 'size' of the respective emission regions to 
be the half hght radius, we calculate the surface density of 
both star formation and molecular gas (which we correct for 
interstellar helium by multiplying our molecular gas masses 
by 1.36). 

Fig. [To] shows our three sources in the context of the 
KS law. The smaller size of the 1.4 GHz emission region 
(respective to the CO emission region) leads to the three 
sources presented in this work lying significantly above the 
~ KS relation, and above the high-z ULIRGS previ- 
ously examined. This indicates an extremely high efficiency 
of star formation respective to the gas density. Even Lock- 
man 38 - which may be forming stars an order of magnitude 
slower than a typical SMG (as per the discussion above) - 
has such an extended CO morphology and compact radio 
emission that its gas surface density is lower than would be 
expected given its SFR, causing it to lie above the high- 
z KS relation. Both HDF132 and HDF254 lie significantly 
(~ 0.8 dex) off the KS slope: it is likely to be a combination 
of the extended CO morphology and compact radio emis- 
sion which conspire to move the sources above the relation. 
As discussed above however, this finding is dependent on 
the assumption of a merger-appropriate CO/H2 conversion 
factor for our sources, which may be significantly underes- 
timating the gas content for our galaxies, particularly the 
disc-like HDF132 which has a gas fraction of 5% assuming 
a ULIRG-appropriate conversion factor. If the smaller con- 
version factor is rejected and instead a Milky Way factor 
is used, then our galaxies fall very close to the high-z KS 



relation of Bouche et al. ( 2007 1 



Using our independent measurements of the MERLIN- 
observed star forming regions and the CO morphology, we 
have also re-calculated Egaa and Esfr for SMGs from the 
literature: these properties are presented in Table|3] For con- 
sistency, all gas surface densities have been adjusted to a 
common CO/H2 conversion factor of 0.8 Mq K km s~^pc^ 
(though an error bar again delineates the possible range of 
values, up to the standard Milky Way value) , and all SMG 
molecular gas masses have had a helium correction factor of 
1.36 applied. These archival SMGs are also shown in Fig. |10| 
There appears to be an offset between the SMGs presented 
in this work (both newly observed and archival) and the KS 
relation for high- 2; star forming galaxies claimed by [Bouchel 
et al. (20071, which suggests that obtaining an estimate of a 



galaxy's gas mass by inverting the canonical KS law can give 
unreliable results. Furthermore, as explained above the off- 
set between our sources and other suggested KS law slopes 
is, perhaps, partly attributable to a combination of uncer- 
tainties in the CO/II2 conversion factor and star formation 
rates for this type of system, but could also provide evidence 
that stars are formed more efficiently in these extreme sys- 
tems. If this is the case, then it may be that a 'universal 
star formation law', which provides a consistent power-law 
relationship between gas content and star formation across 
all systems and environments, is not achievable, or at least 
does not apply in the most extreme cases. 



5 CONCLUSIONS 

We have presented high resolution (sub-arcsecond) detec- 
tions of molecular gas in three z ^ 2 star forming galax- 



ies, two SMGs and one SFRG. We used extracted spectra 
to constrain CO fluxes, and extrapolated total molecular 
gas masses. We analysed the kinematic behaviour of the 
sources revealed in the CO imaging on 0.15" scales, equal 
to ~ 1/3 — 1/2 of the beam FWHM. Our main conclusions 
are as follows: 

• One of our galaxies is unequivocally an early stage ma- 
jor merger (HDF254) with high velocity dispersion and a 
minimal rotational contribution. Lockman 38 appears to be 
a late stage merger, with a low velocity dispersion and a 
smooth morphology, but with significantly disturbed gas 
kinematics revealed by the position-velocity diagram. The 
remaining detected source (HDF132) appears - to the limits 
of our data - consistent with being a disc-like rotating body, 
with double peaked profiles indicative of rotation and a ho- 
mogeneous, non-clumpy morphology. These features may be 
explainable with a late stage merger model, or may be evi- 
dence of a star forming disc. 

• There is a marked discrepancy between the sizes of 
the CO and 1.4 GHz radio emission regions in our high- 
resolution data, with the radio emission region (which traces 
star formation) being more compact than the CO region in 
all cases. This has important implications for our under- 
standing of star formation in high-z starburst galaxies, par- 
ticularly in relation to the role of gas reservoirs, and whether 
starburst galaxies tend to be fuelled by cold fiow accretion, 
or mmor mergers (see [Dekel fc Birnboim|[2006| ) . 

• This size difference leads to an exacerbation of the offset 
of SMGs from the Kennicutt-Schmidt star formation law, 
which has already been shown to be significant (see |Bouche| 
et al.|2007 for a compilation of recent findings in this area). 

• The adoption of a low, ULIRG-like CO/H2 conversion 
factor for SMGs leads to very low gas fractions of < 30%, 
which is far lower than would be expected given their ex- 
treme SFR. The adoption of a Milky Way conversion factor 
(which is somewhat theoretically supported) raises the gas 
mass by a factor of 5. 

• This size difference also manifests in the form of a 
star formation efficiency variation, with 'sub clumps' of our 
sources differing in in efficiency by up to a factor of ~ 5. 
However, the global properties of each source are consistent 
with the canonical L'co-Lfir relation for similarly extreme 
(and sub-mm bright) systems. 

• These are the first PdBI observed 'A' configuration ob- 
servations of high-z ULIRGs for CO(4— >-3), achieving the 
smallest beam size yet for an intermediate-excitation transi- 
tion at a redshift of ~2. The modest exposure times required 
(6 antenna equivalent of ~ 8 hours), the small loss of flux 
(a factor of 1.9 on average), and the detailed morphologies 
and kinematic information revealed suggest that this is a 
profltable approach for similar investigations in the future. 
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